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ABSTRACT 

The 3-D and kinematic structure of the Eskimo nebula, NGC 2392, has been 
notoriously difficult to interpret in detail given its complex morphology, multiple 
kinematic components and its nearly pole-on orientation along the line of sight. 
We present a comprehensive, spatially resolved, high resolution, long-slit spectro- 
scopic mapping of the Eskimo planetary nebula. The data consist of 21 spatially 
resolved, long-sht echelle spectra tightly spaced over the Eskimo and along its 
bipolar jets. This data set allows us to construct a velocity-resolved [Nil] channel 
map of the nebula with a resolution of 10 km s~^ that disentangles the different 
kinematic components of the nebula. The spectroscopic information is combined 
with HST images to construct a detailed three dimensional morpho-kinematic 
model of the Eskimo using the code SHAPE. With this model we demonstrate 
that the Eskimo is a close analog to the Saturn and the Cat's Eye nebulae, but 
rotated 90° to the hne of sight. Furthermore, we show that the main charac- 
teristics of our model apply to the general properties of the group of elliptical 
planetary nebulae with ansae or FLIERS, once the orientation is considered. We 
conclude that these kind of nebulae belongs to a class with a complex common 
evolutionary sequence of events. 

Subject headings: Planetary Nebulae: individual (NGC 2392, NGC 7009, NGC 6543) 
— ISM: kinematics and dynamics — ISM jets and outflows — techniques: spectroscopy 
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1. INTRODUCTION 

Planetary nebulae (PNe) are formed by low to intermediate mass stars. At the end 
of their lives, during the thermally pulsating asymptotic giant branch (AGB) stage, these 
stars expel their chemically enriched outer layers and begin a rapid transformation of their 
degenerate carbon-oxygen core towards their final evolution as a white dwarf. Along this 
route, the previously expelled, cold, neutral, molecular and dusty shell is fully ionized by 
the progenitor star in only several hundred to a few thousand years while the system suffers 
a transition stage known as planetary nebula. Once the central star (CS) leaves the AGB 
the stellar wind gradually develops supersonic velocities and shocks the freshly ionized 
shell, contributing important thermal and hydrodynamic effects to the shell development. 
During this process the nebula also suffers constant changes in terms of optical depth of 
the expanding shell, ionizing conditions and evolutionary shaping takes place. All these 
changes are reflected in the varying spectral emission conditions and morphology of the PN 
as this evolves. 



The Eskimo is a moderately evolved, high excitation PN that presents a roundish 
outline on the sky, although it is well known that this appearance is only due to its polar 
orientation towards the observer since th e real morphology of its m ain shell approximates a 



prolate spheroidal (e.g. Weedman 1968). 



O'Dell fc Balll ( 



19851) and 



Gieseking et al. 



(119851) 



discovered the first bipolar jet in a PN in the Eskimo; a highly collimated, high-speed 
outflow directed nearly towards (and away from) the observer. The bipolar jet is not 
detected in optical images due to its particular orientation and it is only apparent in high 
resolution, spatially resolved spectra. In addition to the former authors, conceptual models 
of the Eskim o based on high resolution spect r oscopy a nd nar r ow-band images hav e been 



presented by 



Reav. 



Phillips fc Cuestal fll999h . and 



Atherton fc Tavloj fll983h. 



Dufour et al. 



Balickl (Il987ch . 



O'Dell et al. 



(1l990h and 



( I2OI2I ). Although with some differences in their 
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models, they all consider a fast expanding, prolate spheroidal inner shell, an outer disk 
and an outer envel ope, the fast bipola r outflow and the knotty structure in the fur of the 



Eskimo. However, 



Zhang et al. 



( 120121 ) have recently argued that NGC 2392 is a twin to 
the tight-waist bipolar planetary Mz-3, but they built their model mostly on arguments 
rather than data, and as we will show here their model does not agree with the most salient 
features present in our data. 



NGC 2392 is also a point and diffuse X-ray source ( IGuerrero et al.ll2005l : iKastner et al. 



2OI2I ). the latter produced by a wind-shock generated hot bubble and likely related to 
the fast expansion of the inner shell. The origin of the X-ray point sourc e is les s clear 
the central star (CS) of the Eskimo has been studied by several authors. 



Heap 



(119771 ) 



classified the CS in the Eskimo as spectral type 06f with Tcs = 35 000 K, a result in 
conflict with that derived from the He II Zanstra temperature that yields ~ 92 000 K. 
This discrepancy lead Heap to suggest that there must be a second hotter star responsible 
for the nebular ionization. Several autho rs have calculated an effective temperature , 



Tpff, in the range 40 000 K - 45 000 K. (iMendez et al. 



Kudritzki et al. 



2011 



Pauldrach et al. 



2004 



19971 ). However this temperature is not high enough to explain the high 



stages of ioniz ation from some ions such as O IV and Ne V which have been o bserved 



in the nebula ( iPottasch.Surendiranath &: Bernard-Salas 



2008 



Natta et al 



19801). Other 



studies have suggested that the effect i ye temperature of the CS must be around 74 00 



to 80 000 K (iTinkler fc Lamers 



2002 



Pottasch.Surendiranath &: Bernard-Salas 



20081 ). 



Pottasch.Surendiranath fc Bernard-SalasI (|2008h fi nd that the central star must have evolved 



from a 1.7 M0 progenitor. iCiarduUo et al. 



(119991 ) detected a possible fa int companion t o 



the C S with HST in the I band but invisible in the V band. Recently, 



Danehkar et al. 



( I2OI2I ) used photoionization models to estimate that the companion star must have a 
Teff = 250 000 K, which is much higher than any other value proposed earlier. Although 
the presence of a faint and hot companion to the CS of the Eskimo is likely, there is no firm 
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or unambiguous evidence of its existence to date. 



The most recent distance estimates to t 



2008h with d = 1.6 kpc±130 pc. 



le Es kimo nebula are those of lStanghelhni &: Shaw 



Pottasch. Surendiranath fc Bernard-Salad ( 1201 ll ) with d = 1.8 kpc. 



Frewl (2008) with d = 1.28 kpc ±130 pc and 



The spectroscopic observations reported here were obtained with the aim to carefully 
disentangle the different kinematic components of the Eskimo nebula and their relative 
location in space . We use the morpho-kinematical modeling code SHAPE to construct a 
detailed three dimensional model of the Eskimo from the kinematic data combined with 
high spatial resolution images from HST. A 3D model is produced as the final output that 
clearly illustrates the structure of the Eskimo nebula. In addition, we identify significant 
structural similarities between the Eskimo (NGC 2392), the Saturn (NGC 7009) and the 
Cat's Eye (NGC 6543) nebulae, and in general with the family of elliptical PNe with ansae 
or FLIERS, that provide a fresh unifying view of this group of PNe. In §§ 2—3 we first 
describe the observations, data reduction, and the various morpho-kinematic structures 
of the nebula. Isovelocity channel maps are presented in §4. Then we describe the 3-D 
SHAPE model of NGC 2392 and a comparison of the Eskimo with the Saturn and Cat's 
Eye nebula (§5). Finally in §6 we summarize the conclusions of this study. 



2. OBSERVATIONS AND DATA REDUCTION 



Long-slit, echelle, spectroscopic observations of NGC 2392 were obtained at the 
Observatorio Astronomico Nacional at San Pedro Martir, (SPM), Baja Califor nia, Mexico, 



20031 ) on the 2.1 



with the Manchester Echelle Spectrometer (MES-SPM) ( iMeaburn et al.l 
m telescope in its //7.5 configuration. The observing run took place in 2002, January 7 - 
10. MES-SPM was equipped with a SITE-3 CCD detector with 1024 x 1024 square pixels. 
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each 24 /im on a side (= 0.312 arcsec pixel"^). We used a 90 A bandwidth filter to isolate 
the 87th order containing the Ha and [N II] nebular emission lines. Two pixel binning was 
employed in both the spatial and spectral directions. Consequently, 512 increments, each 
0'.'624 long gave a projected slit length of 5^32 on the sky, except for slit position g where 
we used no binning (i.e. 1024 x 1024). For the majority of the exposures, we used a 70 fim 
(=0.95") slit, giving a velocity resolution of 9.2 km s^^ (= 0.312 arcsec pixel"^ and 5.6 
km s~^ for slit g) and for two exposures, slits a & r, we used 150 yum (= 1'.'9 and = 11.5 
km s~^). We obtained 18 consecutive and tightly spaced positions over the Eskimo with 
P.A. = 0°, one position with P.A. = 110° across the center of the nebula, slit s, and two 
positions with P.A. = 70°, slits t & u. Slit u sits on the bipolar jet. The slit positions are 
indicated and labeled in Figure 1 on a WFPC 2 image of the Eskimo obtained from the 
HST archive. All spectra were acquired using exposure times of 1800 s. In order to establish 
the exact position of the slit in each pointing, the slit position on the sky was recorded with 
an automatic procedure available in MES-SPM prior to the spectroscopic exposure. 

The data was reduced by using standard IRAfQ tasks to correct bias and remove 
cosmic rays. The spectra were calibrated in wavelength against the spectrum of a Th/Ar 
lamp to an accuracy of ±1 km s~^ when converted to radial velocity. All spectra presented 
in this paper are corrected to heliocentric velocity (Vhei)- 



All the spectroscopic data have been drawn 



matic Catalogue of Galactic Planetary Nebulae" flLopez et al 



rom "The San Pedro Mdrtir Kine- 



20121 ) and available at 



http://kincatpn.astrosen.unam.mx In the last part of this work we also make use of spectra 



for the Saturn (NGC 7009) and Cat's Eye (NGC 6543) nebulae, these have also been drawn 



^IRAF is distributed by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science foundation 
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from the same source, they both have been observed with the same detector as described 
above, binned 2x2 and with a 150 yum sht (= 1'.'9 and = 11.5 km s~^). The spectra for the 
Saturn and Cat's Eye nebulae are shown in Figure 9. 

Additional steps where taken to produce velocity cubes in [N II] 6584 . These were 
implemented using FORTRAN routines. The additional tasks were: a) Continuum emission 
was removed by fitting a quadratic function to each row of each two-dimensional spectrum, 
b) An astrometric solution was employed for each of the image+slit using several stars near 
the Eskimo in order to find the exact slit position of each exposure. The position of each 
star was taken from Sky View data, c) The data were photometrically calibrated using the 
[N II] 6584 HST image. In this way, the observed slit spectra oriented N-S were combined 
and interpolated to const ructed isovelocity chan n el ma ps of the Eskimo. This techniques 



are described in detail in 



Garcia-Diaz &: Henneyi (120071 ). 



The resulting data cubes for the [N II] 6584 line are shown in Figure 5 as isovelocity 



Balick et al. 



(jl987al ) for a similar 



channel maps, each isovelocity map is 60 km s~^ wide (see 
data representation with narrower velocity intervals.). The derived moment maps for the 
same emission line are shown in Figure 6. 



3. KINEMATICS 

The bi- dimensional [N II] emission line spectra or position-velocity (P-V) arrays for 
all individual slit position are shown in Figures 2 and 3 (except for slit position r since the 
[N II] line emission was extremely faint at this location). For each slit position we show a 
couple of P-F arrays: the observed [N II] P-F array is on the left panel and on the right 
panel is the corresponding synthetic P-F array derived from the SHAPE model. Spatial 
offsets are in arcseconds. The stellar continuum from the CS has not been subtracted and 
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this is apparent in the P-F arrays corresponding to slit positions j, s and u. All spectra 
are presented within a heliocentric velocity range of —135 to +278 km s~^. We choose this 
range because it allows to appreciate the full complex structure of the line profiles and 
clearly identifies the beginning and the end of the high-velocity, bipolar, coUimated outflow. 
We derive a heliocentric systemic velocity, V^ys= 70.5 km s~^ by using slit position j, which 
passes through the central star. 

In order to facilitate the discussion of the kinematic information, we have labeled in 
Figure 4 each emission region seen in the P-V arrays in slit positions j and u. Each of these 
components are present throughout the spectra presented in Figures 2 and 3 and are the 
basis of the interpretation and our model of the Eskimo. A description of each component 
is given below, please refer to Figure 4 for the remainder of this section. 

3.1. THE INNER SHELL 

The inner shell is recorded from slit positions g to m that cut the nebula at P.A. = 0° 
over the central region and also in slits u and s (P.A. = 70°). The line profiles from the inner 
shell corresponds to a tilted and distorted (peanut-like) velocity ellipse ( VE) with bright 
regions over its outline that are part of the filamentary structure projected over the surface 
of this shell. In the central slit, j, the velocity ellipse covers a velocity range Vhei — —50 to 
190 km s~^, considering that Vsys— 70.5 km s~^, then the expansion velocity of the inner 
shell along the line of sight is Kxp— ±120 km s~^, this is probably the largest expansion 
velocity for an elliptical shell in a PN and presumably corresponds to the expansion of the 
prolate ellipsoid, indicated by the shape of the P-F array, close to its major axis. Note that 
the tilt of the VE is more pronounced in sht g than in sht m, that correspond to the edges 
of the shell. There is also a small velocity shift in the limits of the VE's corresponding to 
the contiguous positions on either side of the central slit. This shift cannot be read out 
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from the P-F arrays in Figure 2 due to the scale, this shift is of the order of 10 km s~^; 
for sht i the VE is shghtly redshifted and for sht k the shift is in the opposite sense. These 
elements indicate that the main axis of the fast-expanding ellipsoidal shell is pointing not 



pole on) but slight ly south and to the southwest. This is in 
O'Dell et all (|l99o|) that derived a P. A. = 202° for the main 
axis of the ellipsoid or prolate spheroidal shell and an inclination ~ —19° with respect to 



directly towards the observer 
accordance with the model of 



the line of sight. In our SHAPE model (see below, section 5) we obtain P. A. = 205° and 
an inclination ~ —9° with respect to the line of sight. 



3.2. THE OUTER SHELL AND COMETARY KNOTS 

The outer shell is revealed throughout all the P-F arrays in Figures 2 and 3 as the 
central and narrow velocity ellipse expanding with a velocity \4xp — 16 km s~^. This 
shell is usually referred to as the hood of the Eskimo. The shape of the VE's is smooth 
and uniform, indicating that the outer shell is close to an expanding sphere. Some of the 
cometary knots that are seen in projection against the hood of the Eskimo are apparent as 
bright knots on the extremes of the narrow VE of the outer shell throughout the different 
positions where the slits intersects the knots or their bright tails, and their velocities 
coincide with V^ys, which means that either these knots are inert or any motion in the 
system of these knots must be tangential to the sightline and implying that they must be 
distributed in a disk-like configuration. Since the major axis of the inner ellipsoidal shell is 
nearly pole-on that means that the disk must be close to the plane of the sky and likely 
coinciding with the equatorial region of inner shell. 
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3.3. THE CAPS 



In addition to the inert cometary knots there is an extended system of bright knots and 
diffuse material located in the outer shell that make-up the fur in the Eskimo's hood. These 
knots are located below and above the VE that describes the inner shell and shifted in 
velocity with respect to the thin VE that represents the outer shell. We call these systems 
or groups of knots simply caps. Their projected location corresponds to the fur of the 
Eskimo's hood, as the cometary knots described above, but their radial velocities indicate 
that they are distributed in three extended groups. The knots located in the northern 
section all show only receding velocities whereas the group in the southern section show 
both receding and approaching radial velocities but with a different spatial distribution. 
Those with receding radial velocities appear at projected locations slightly above the knots 
with blueshifted radial velocities. Considering that the tilt for the inner shell described 
above applies also to the rest of the fur structure, this projected spatial segregation of the 
receding and approaching system of bright knots in the southern section indicates that 
we are looking at caps in front and behind the inner shell and with an expansion velocity 
~ ±55 km s^^. In the northern section there is only the receding cap and presumably 
located in the back part of the nebula. The location and approximate extent of the caps 
is indicated in Figure 4. These caps very much resemble the system of FLIERS and 
low ionization emission regions, like those in NGC 7662, NGC 7009 and NGC 6543 (e.g 



Balick et al 



1987b 



1998 



Reed et al. 



19991), see Figure 9. 



3.4. THE COLLIMATED, HIGH VELOCITY, BIPOLAR OUTFLOWS 

The bipolar jets are clearly appreciated in slit u. They are seen to be launched right 
from the star (the bipolar jets emerge from the stellar continuum) at ~ Vhoi — 95 and 
+235 km s~^; considering that Vgys = 70.5 km s~^, they are launched from the source at 



- 11 - 



V^xp = ±165 km s~^. The jets continue increasing their velocity with distance to reach 
top velocities Vhei — HO and +250 km s~^or Kxp— ± 180 km s~^. At the point where the 
bipolar jets seem to emerge from the inner shell, approximately 12 arcseconds away from 
the stellar continuum, the jet emission suffers a discontinuity and its full extent gets split in 
two sections, in both directions, at this point both jets seem to suffer a slight acceleration. 
All jet components have similar brightness as well as spatial and velocity extents. The 
discontinuity in the bipolar jets may arise from a sudden change of environment, when 
moving from the inner hot shell into the realm of the outer shell or to an episodic 
ejection. The presence of the bipolar jets is also apparent from shts c - q at the locations 
where the slits intersect the jet. The redshifted fast velocity components are observed from 
slits c to i and the blue-shifted counterparts from shts k to q. The projected width of the 
jets in these slits is of the order of 10 arcseconds and remains remarkably constant along 
the entire length of the jet. 

4. ISOVELOCITY CHANNEL MAPS 

In order to provide an alternative overview of the kinematic structure of the Eskimo we 
constructed isovelocity channel maps of the [N II] optical emission hne. The channel map 
representations are very useful to reveal large-scale, spatially coherent features that are not 
readily apparent in the individual long-slit spectra. The maps were derived from slits a - r, 
each velocity cube is an image of the nebula within a 60 km s~^ range in velocity, extending 
from -1-240 to —110 km s~^. These isovelocity channel maps are presented in Figure 5, 
where relative positions are given in arcseconds with respect to the CS. The velocity interval 
and limits have been chosen to produce a clear snapshot of the main kinematic structures 
of NGC 2392. The first (a) and last (e) panels in the figure clearly show the coUimated 
bipolar outfiow at the extreme velocities covered by the isovelocity channel maps. The 
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second (b) and fourth (d) panels reveal the location of the redshifted and blueshifted caps 
outside the inner shell, note the slight but apparent spatial displacement of the inner shell 
between the redshifted (b panel ) and the blueshifted (d panel) maps due to its tilt with 
respect to the line of sight. Finally, the central (c) panel shows the structure of the Eskimo 
at velocities close to systemic; cometary knots and caps are clearly discerned comparing the 
central panels. 

We have also constructed the kinematic moment maps for the same [N II] emission 
line. These maps of derived parameters offer a complementary approach to the isovelocity 
channel maps. The parameters shown in Figure 6 are the total line surface brightness, 
SB, mean heliocentric velocity, Vhei, root-mean-square (RMS) velocity, a, and skewness Sk- 
These quantities are defined in terms of the velocity moments, Mk of the hne profiles I{v): 

Mfc= / v^I{v)dv (1) 

and the limits of the integration are vi = —140 km s~^ and V2 = +280 km s~^, which 
represent a wide enough range to include most of the emission of interest. Thus, the zero 
moment, Mq, 

SB^Mo^ / I{v)dv (2) 
The first velocity moment. Mi is, 

rv2 

Ml = / vl{v)dv, (3) 

J VI 

and the mean heliocentric velocity is calculated then from the zero and first moment. 



The root-mean-square velocity width, a, is derived from the second moment, M2, 



V2 

M2= I v^I{v)dv 



(5) 



then. 



FinaUy, from the third moment the velocity skewness map is derived: 



M- 



V2 

3 = / v^I{v)dv 

V\ 



(7) 



and. 



Sk = M3/M0 (8) 

The first panel in Figure 6 shows the map of the [N II] line intensity distribution or 
spatial variations over the nebula within the full velocity range of integration. The second 
panel displays the distribution of Vhei that highlights the regions corresponding to the 
shells and caps. The RMS velocity is displayed in the third panel, since a is sensitives to 
low emission gas at high velocities , the bipolar coUimated outflows become apparent here, 
finally the last panel contains the skewness map that shows that the velocity field within 
the nebula is mostly symmetric. 
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MORPHO-KINEMATIC MODELING WITH SHAPE 



The 3-D morpho-kinematic structure in th e [N II] 6584 emis s ion o. 



mode 



ed using the code SHAPE, developed by 



Steffen fc LopezI fl2006f ) and 



NGC 2392 has been 



Steffen et al 



( 2010[). Similar analyses have been pe rformed recently on the planetar y nebulae NGC 63 37 



Garcia-Diaz et al. 



and Hb 5 (ILopez et al 



2009h. NG C 7009 flSteffen et al. 



2012ah . 



20091 ). NGC 6751 ficiark et al.ll2010h . 



Modeling with SHAPE follows three main steps. First, defining the geometrical forms 
to use; SHAPE has a variety of objects such as a sphere, torus, cone, cube, etc. whose 
basic forms can be modified by the user (e.g. squeeze, twist, boolean, etc). Second, an 
emissivity distribution is assigned to each individual object or structure, and third, a 
velocity law is chosen as a function of position. SHAPE gives as result a two dimensional 
image and synthetic P- V arrays, which are rendered from the 3D model to be compared 
with the observed data. The parameters of the model are then iteratively adjusted until a 
satisfactory solution is obtained. 

The Eskimo has been modeled using the [N II] i/ST image and the P-F spectra shown 
in Figures 2 and 3. Our model includes five main features: the inner shell, the outer shell, 
the cometary knots, the caps and the high-velocity bipolar outflows. We modeled the 
outer and inner shells starting with a sphere and then modified them to conform with the 
observed image and P-F arrays. For the outer shell we made no effort to deform a sphere 
since there is no information on the data to guide us in this task, its real shape may slightly 
differ from a sphere but this is of no consequence for the global model and for this reason 
we call it an spheroid instead of only a sphere. For the inner shell, strings of filaments were 
also added in its inner sections and over its surface (see Figure 7). 

For the velocity field we started with a homologous expansion, i.e. radial velocity 
vectors linearly increasing in magnitude with distance and added a poloidal component. 
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as described in ISteffen et al.l (120091 ) and iLopez et al.l (l2012al ). however it was found that 



a velocity field with only radial velocity components is sufficient to describe the velocity 
law in this case. The stationary cometary knots where modeled individually from cylinders 
with emissivity variations along the axis of the cylinder and distributed over a ffattened 
toroid or disk sitting near the equatorial region of the inner shell. For the filamentary and 
diffuse large scale emission pertaining to the caps we used again sections of spheres with 
emissivity variations along their axes distributed over sections of a fiattened torus placed 
at the border of the outer shell and projected in front and behind the tilted inner shell. 
The bipolar coUimated outfiows were built from long cylinders whose length was arbitrarily 
limited to the size covered by the set of spectra. 

The results of this process are the synthetic P- V arrays shown next to the observed 
ones in Figures 2 and 3. The main uncertainty in the model is the size and the placement 
of some structures along the line of sight since the information along the third dimension or 
into the plane of the sky is derived from the spatially resolved spectra assuming that the 
nebula is expanding in a homologous way. For example, the bipolar jets may have a larger 
extent than observed or the caps may be located either at the border of the outer shell, as 
in NGC 6543, in between the inner and the outer shells or even closer to the inner shell, as 
in NGC 7009, (see Figure 9). For these reasons our global solution is not unique in these 
details, however, since in this case there is a very good spectral coverage and the nebula 
has so many independent complex elements, the solution is very well constrained by the 
excellent match of the synthetic P-F arrays to the actual data. 

The resultant SHAPE mesh model for the Eskimo, rotated 90°to the line of sight, 
before rendering is shown in Figure 7, where the individual main components are labeled. 

The results of the final rendered model are shown in Figure 8, where the top panels 
show first the synthetic P-F array for the central position, slit j, then the synthetic image 
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from the model for the Eskimo as seen on the sky and next to it the same model image 
but rotated 100° into the plane of the sky and also rotated clockwise 45°. A composite 
HST image of the Cat's Eye nebula, NGC 6543 is next presented for comparison. In the 
bottom panels we have used the SHAPE model from Steffen et al. (2009) for NGC 7009, 
the Saturn nebula. That model has been rotated 90° for this work to show how its central 
line profile would appear if we were looking at it pole-on, as in the case of the Eskimo. The 
corresponding model image is next and to its right is the model image of NGC 7009 as seen 
on the sky. Finally an RGB HST image of the Saturn nebula is shown in the last panel for 
comparison. The structural similarities between the Eskimo, the Cat's Eye and the Saturn 
nebulae are apparent and striking. 

As an additional point of comparison among these PNe, in Figure 9 are presented the 
observed P-F arrays for NGC 6543 and NGC 7009 obtained along their major axes, the 
data have been dr awn from the SPM kinematic catalogue of galactic planetary nebulae 



( iLopez et al. 



20121 ). together with a synthetic P-F array of the Eskimo obtained also 
along its major axis after rotating the model nebula 90° to derive the corresponding 
bi-dimensional line profile. The velocity scale for the synthetic P-F array has been squeezed 
by a factor of two and blurred with a gaussian filter to make it comparable with the rest. 
It is clear that all these P-F arrays (nebulae) share similar structural elements. 



6. DISCUSSION AND CONCLUSION 

Our aim in this work has been to disentangle in detail the main morphological and 
kinematic elements of the complex Eskimo nebula, NGC 2392. In order to present the 
results in an easy to visualize form we have paid particular attention to relate key structures 
in the projected morphology of the Eskimo from HST images with spatially resolved, 
long-slit, echelle spectra. Channel maps and kinematic moment maps have been created 
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from the set of long-slit spectra to aid visualize the kinematic structure of the Eskimo. This 
relation between morphology and kinematics has been modeled using the code SHAPE in 
order to produce a 3D model of the object. We find that many of our results are in general 
agreement with early works on the Eskimo, cited in the introduction, that gave the first 
insights into the actual st ructure of th i s PN. However, our data and model do not agree with 



the recent suggestion by iZhang et al. 



(120121 ) that consider the Eskimo twin of Mz-3, 
a bipolar nebula with a very tight waist and a complex outflow structure external to the 
main lobes. The present work presents the largest published collection of spatially-resolved, 
long-slit, echelle spectra for the Eskimo that together with the SHAPE model provide a 
first, clear, data-based 3D representation of this complex pole-on oriented nebula. 

The data and modeling presented here yields the following results for the different main 
components of the Eskimo: the outer shell is an oblate spheroid with Vc^p = 16 km s^^. 
This is probably the only structural element that preserves some memory of continuous 
evolution from the end of the AGB stage up to now. For an outer shell angular radius r = 
23", Vexp= 16 km s~^ and adopting a distance D = 1.4 kpc, the estimated mean kinematic 
age would be 9300 years. 

The inner shell is a deformed (peanut-like) prolate spheroid with V^xp = 120 km s~^ 
along its major axis. Its orientation is nearly pole-on, tilted 9° with respect to the line 
of sight and pointing towards P. A. = 205°, i.e. to the south and west. Multiple bright 
filaments are distributed over its surface that under ground-based seeing conditions give 
it the appearance of a continuous bright border rim. Its side view shape from the model 
(see Figures 8 and 9) is very similar to the shape of the inner shells in NGC 6543 and NGC 
7009. The velocity law in our model assumes a homologous expansion with a Hubble- type 
velocity law of the form v = k ■ r /vq, where is a constant, r is the distance from the center, 
and tq is the distance at which the velocity k is reached. The model does not require a 
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poloidal velocity component in order to match its current expansion conditions and shape, 
only different values of k are applied to the different components of the nebula . If the 
outer and inner shell formed at approximately the same time, at the current expansion 
rate of the latter along its major axis it would have overtaken the outer shell in only 
1250 years. Therefore, either the inner shell had a much slower acceleration in the past 
or its present expansion rate has been in operation for a relatively short time compared 
to the development of the outer shell. Although the current exact separation between the 
approaching tip of the inner shell and the border of the outer shell cannot be discerned, our 
model assumes that the inner shell has not reached the border of the outer shell yet. The 
fast expansion of the inner shell along a preferred symmetry axis occurred at a much later 
time than the exit from the AGB stage, once the envelope was fully ionized and the central 
star had developed a sufficiently fast wind to influence the expansion of the shell through 
the mechanical and thermal energy produced by the shocked wind on the previously 
expelled matter. Once the fast expansion of the inner shell started, this clearly took a 
non -isotropic mode . The inner shell harbors a hot bubble with extended X-ray emission 



but 



Guerrero et al. 



(|2005[ ) point out that the stellar wind velocity does not seem to be high 
enough to account for the X-ray luminosity and briefly discuss the possible contributions 
from a binary companion and/or the fast bipolar outflow to resolve this discrepancy. 



The bipolar jet is detected emerging right from the star at V^^p = 165 km s^^, it has 
a near pole-on orientation, with its axis tilted some 5° further to the southwest than the 
major axis of the inner shell; the jets get split into two sections as they leave the inner shell 
and then reach Vexp = 180 km s^^, the jets show an approximately constant width all along 
its detected length, indicating a high degree of collimation. Its projected length is 25" but 
tilted some 30°with respect to the line of sight, assuming a constant velocity, the jets have 
been expanding for approximately 1800 years, without considering possible episodic events. 
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The caps are distributed in three main extended sections that in projection form 
most of the fur of the Eskimo's hood. Each of these sections has a bright portion that 
is conspicuous in the fur of the Eskimo's hood, but the real extent of each of them is 
larger than the bright portions as indicated by the spectroscopic information. Two of these 
sections are in the back part, mainly to the north and west of the inner shell and show the 
corresponding receding velocities, the third section is in the front and to the south of the 
inner shell and shows the expected approachi ng velocity: the caps are expanding at 



55 km s ^ . Thes e caps correspond to what 



O'Dell et al. 



O'Dell k Balll (119851 ) called wisps and later 



( I2OO2I ) called the north, west and south fuzz. The caps are not aligned with the 
axis of the bipolar jet, a similar situation occurs in NGC 6543 (see Figure 9). Considering 
that the inner shell expands twice as fast as the caps and these seem to be located beyond 
the tips of the inner shell, this suggests that either the caps (perhaps together with the 
equatorial disk of cometary knots) were part of a pulse of mass loss that occurred sometime 
after the formation of the PN or are part of the outer shell. In our model the caps are placed 
at the edge of the outer shell, if their expansion velocity has been constant through time, 
it would take them 2700 years to reach the outer envelope from the star and, as mentioned 
before, the inner shell would have caught up with them in half that time. The reasonable 
alternative is then that they are originally part of the outer envelope. The structure of the 
caps in our model is very similar to the caps in NGC 6543 (see Figure 9). 

The cometary knots are distributed over an equatorial disk, external to the inner shell 
at a radius of ~ 17" that projects them over the fur of the Eskimo's hood and mixes them 
with the caps; these knots are at the systemic velocity, i.e. they seem inert or with no 
radial velocity component, it is likely that they are slowly expanding radially with respect 
to the central star and have bee n formed clo s e to i ts current location, or near the main 



ionization front, as discussed by 



O'Dell et al 



( 12002[ ). Note that a similar equatorial knotty 



ring is present in NGC 6543 (see Figure 9) and traces of a similar structure are discernible 
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in NGC 7009. 

Our model has lead us in a natural way to compare the striking similarities of the 
Eskimo nebula with NGC 6543 (the Cat's Eye) and NGC 7009 (the Saturn) nebulae, and 
by extension with the rest of the group of elliptical n e bulae with ansae or FLIERS (e.g. 



NGC 3242, NGC 6826 and NGC 7662). 



Balick et al. 



( 1l998l ) have described in detail the 



characteristics of these PNe, they all show a closed inner shell with a bright rim, an outer 
shell with low ionization emission regions and a coUimated bipolar outflow. Comparing the 
images in Figure 8 it is straightforward to relate the curious, peanut-like, shape of the inner 
shell in the Eskimo with the inner shells in NGC 6543 and NGC 7009. An equatorial disk of 
cometary knots as the one described for the Eskimo is apparent in the image of NGC 6543 
and the structures related to the caps in the Eskimo can also be clearly identified in NGC 
6543 and NGC 7009, though in the latter they seem to be distributed over a cylindrical 
shell right outside the inner shell. In both cases they also have bipolar collimated outflows. 
With these analogies in mind it becomes easy to visualize how the Cat's Eye nebula and the 
Saturn nebula would look like if viewed pole-on, very much like the Eskimo. 



Steffen et al. 



(120091 ) have made a model of the Saturn nebula and its pole-on view is shown in the second 



frame of the lower panel of Figure 9. 

This group of ellip tical nebulae with ansae or FLIERS is particularly interesting since 



as 



Kastner et al. 



( 120121 ) show, all the members of this group are among the relative small 
number of PNe where extended X-ray emission has been detected. In addition, they all 
lie at distances from the galactic plane that range from ^ 450 - 800 pc, exceeding the 
PN scale-height of 250 pc for elliptical PN. This suggests that the progenitors are old-disk 
stars of relatively low mass but still massive enough so their post-AGB wind evolve at an 
adequate pace to be able to develop a closed inner shell and a hot bubble. If the remnant 
stellar mass is too low, the stellar wind will accelerate at a very slow rate, keeping the 
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initial two-wind interaction in a momentum conserving state and the shocked wind region 
remains isothermal, allowing the shell to evolve and expand, thus never reaching conditions 
to form a hot. X-ray emitting bubble. 

Although no binary companion has been found in the central star of the Eskimo, 



strong argumen 


ts in 


(Guerrero et al. 


2005 



Kastner et al 



2OI2I ) cannot be explained by the the 



relatively low stellar wind velocity, Voo — 400 km s ^ and low effective temperature 



-efl 



45 X 10'^ K of the central star; the star has a mass loss rate M < 3.2 x 10'' 



Mq yr ^, (IKaschinski. Pauldrach &: Hoffmann 



2OI2I ). The characteristics of the central 



star of NGC 2392 are also inconsistent with the high stages of ionization observed in the 
nebula or the signific ant Zanstra discrepancy derived from the photoionization models (e.g. 



Danehkar et al. 



2OI2I) 



Furthermore, the high expansion velocity of the inner shell seems to require a driving 
force beyond the modest characteristics of the central star since it is only the slightly higher 
thermal press ure of the hot X-ray emitting gas over that of the inner shell that drives 



its expansion ( iGuerrero et al. 



20051 ). and there is still the fast bipolar collimated outflow 
that needs a plausible launching engine. All these elements favor the presence of a binary 
companion to the central star of NGC 2392. 

Just as in the case of the Eskimo, there is no firm evidence of binary n uclei in the oth er 
memb ers of the group of elliptical nebulae with ansae or fliers. Recently, 



Douchin et al. 



( 2012 ) list NGC 6826 as a possible binary with a period of 0.619 days, but 



Jevtic et al. 



( l2012[ ) after a year of continuos observing with Kepler find only stochastic brightness 



variations in the central star of NGC 6826. 
these brightness variations found by 



Prima, Massa fc Cantielld (j2012[ ) consider that 



Jevtic et al. 



mm in the central star of NGC 6828 



are possibly related to subsurface convective layers, just as the ones they have found in the 
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central star of NGC 6543 and where they derive a high rotational velocity for the central 
star. High stellar rotation in the central star of a PN may be an indication of spin-up by a 
past merger, tidal interaction or mass accretion in a binary system. High rotational stellar 
velocities help also in the pro duction of toroidal magnetic fields tha t serve as launching 



platforms for bipolar jets (e.g. 



Garcia-Segura. Lopez &: Franco 



20051) 



A detailed reexamination of the evolutionary history, the possible common mechanisms 
and sequence of physical events in the formation of the members of this group that make 
them share so many outstanding characteristics is evidently called for, but out of the scope 
of the present paper. 

We have presented a detailed morpho-kinematic analysis of the Eskimo planetary 
nebula, NGC 2392. Our resultant 3D model reveals that this nebula is a close analog to the 
Cat's Eye (NGC 6543) and the Saturn (NGC 7009) nebulae once the orientation effect is 
considered. 
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Fig. 1. — Location of each slit position is indicated and labeled on an HST [N II] 6584 image 
of the Eskimo. North is up, east left. 
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Fig. 2. — Mosaic of bi-dimensional P-F arrays labeled according to slit position. For each 
slit position we show a couple of P-F arrays: the observed [N II] 6584 P-F array is on the 
left panel and the corresponding synthetic P-F array derived from the model is on the right 
panel. This figure shows slit positions a-j, the remainder positions are shown in Figure 3. 
The top of the array corresponds to the north end of the slit. 
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Fig. 3. — As in Figure 2 but for the slit positions k - u. The top of the array corresponds 
to the north end of the sht for sht positions k - q. For sht position s the top of the array 
corresponds to the northwest end of the sht and for sht position t and u the top of the array 
to the northeast end of the sht. 




Fig. 4. — Finding chart for individual emission regions seen in the F-F arrays in sHt positions 
j and u. Each of these components is present throughout the spectra presented in Figures 2 
and 3 
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Fig. 5. — [N II] 6584 emission- line channel maps constructed from the long-slit echelle 
troscopy. Each isovelocity image is 60 km s~^ wide. 
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Fig. 6. — Moment maps for the [N II] 6584 line. Panel a) Surface brigthness; b) mean 
heliocentric velocity; c) rms velocity and d) skewness. All moment maps are integrated over 
the full velocity range, VJiei^— 140 to +280 km s~^, but the scale bars at the top of the 
maps show only the interval for the most representative values displayed within the nebula, 
excluding the background 
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Fig. 7. — SHAPE mesh model of the Eskimo, rotated 90° to the hne of sight, before rendering. 
The individual main components are labeled. 
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Fig. 8. — Upper panels, from left, The synthetic line profile for the central slit j. The 
Eskimo image, as seen on the sky, derived from the model. The previous image rotated 100°, 
into the plane of the sky and 45° clockwise. A composite Chandra + HST image of NGC 
6543, shown for comparison. Lower panels: A synthetic line profile from a slit located at 
the center of a pole-on view of NGC 7009, from the SHAPE model of Steffen (2009). The 
corresponding model image of NGC 7009 as it would appear if seen pole-on. The model 
image of NGC 7009 as seen on the sky. HST image of NGC 7009 shown for comparison. 
Image credits for the HST images: NGC 6543, J. P. Harrington & K. J. Borkowski. NGC 
7009, B. Balick, NASA/HST. 
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Fig. 9. — Left and central panels, the observed line profiles for a central slit in NGC 6543 
and NGC 7009, respectively. Right panel, the synthetic line profile for NGC 2392 showing 
how it would like along the axis from a side view, as the former cases. 



